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Abstract 

This is the first in a series of papers on the construction of explicit solutions 
q | to the stationary axisymmetric Einstein equations which describe counter-rotating 

disks of dust. These disks can serve as models for certain galaxies and accretion disks 
in astrophysics. We review the Newtonian theory for disks using Riemann-Hilbert 
methods which can be extended to some extent to the relativistic case where they 
lead to modular functions on Riemann surfaces. In the case of compact surfaces 
^ . these are Korotkin's finite gap solutions which we will discuss in this paper. On 

the axis we establish for general genus relations between the metric functions and 
hence the multipoles which are enforced by the underlying hyperelliptic Riemann 
surface. Generalizing these results to the whole spacetime we are able in principle 
to study the classes of boundary value problems which can be solved on a given 
Riemann surface. We investigate the cases of genus 1 and 2 of the Riemann surface 
in detail and construct the explicit solution for a family of disks with constant 
angular velocity and constant relative energy density which was announced in a 
previous Physical Review Letter. 

PACS numbers: O4.20.Jb, 02.10.Rn, 02.30.Jr 

1 Introduction 

The importance of stationary axisymmetric spacetimes arises from the fact that they can 
describe stars and galaxies in thermodynamical equilibrium (see e.g. §])• However the 



1 



complicated structure of the Einstein equations in the matter region which are apparently 
not completely integrable has made a general treatment of these equations impossible up 
to now. Thus only special, possibly unphysical solutions like the one of Wahlquist were 
found (in M it was shown that the Wahlquist solution cannot be the interior solution for 
a slowly rotating star). Since the vacuum equations in the form of Ernst are known 
to be completely integrable || |7|, the study of two-dimensional matter models can 
lead to global solutions of the Einstein equations which hold both in the matter and in 
the vacuum region: the equations in the matter, which is in general approximated as an 
ideal fluid, reduce to ordinary non-linear differential equations because one of the spatial 
dimensions is suppressed. The matter thus leads to boundary values for the vacuum 
equations. 

Disks of pressureless matter, so-called dust, are studied in astrophysics as models for 
certain galaxies and for accretion disks. We will therefore discuss dust disks in more 
detail, but the used techniques can in principle be extended to more general cases. In the 
context of galaxy models, relativistic effects only play an important role in the presence of 
black-holes since the latter are genuinely relativistic objects. A complete understanding 
of the black-hole disk system even in non-active galaxies is therefore merely possible in 
a relativistic setting. The precondition to construct exact solutions for stationary black- 
hole disk systems is the ability to treat relativistic disks explicitly. In this article we will 
focus on disks of pressureless matter. By constructing explicit solutions, we hope to get a 
better understanding of the mathematical structure of the field equations and the physics 
of rapidly rotating relativistic bodies since dust disks can be viewed as a limiting case 
for extended matter sources. Hence we will discuss relativistic effects for models whose 
Newtonian limit is of astrophysical importance. We will investigate disks with counter- 
rotating dust streams which are discussed as models for certain SO and Sa galaxies (see || 
and references given therein and [ IP| , [TT|]). These galaxies show counter- rotating matter 
components and are believed to be the consequence of the merger of galaxies. Recent 
investigations have shown that there is a large number of galaxies (see ||, the first was 
NGC 4550 in Virgo) which show counter-rotating streams in the disk with up to 50 % 
counter-rot at ion . 

Exact solutions describing relativistic disks are also of interest in the context of numerics. 
They can be used to test existing codes for stationary axisymmetric stars as in |12. 13| . 



Since Newtonian dust disks are known to be unstable against fragmentation and since 



numerical investigations (see e.g. [|T4|) indicate that the same holds in the relativistic case, 
such solutions could be taken as exact initial data for numerical collapse calculations: due 
to the inevitable numerical error such an unstable object will collapse if used as initial 
data. 



In the Newtonian case, dust disks can be treated in full generality (see e.g. (T^j) since 
the disks lead to boundary value problems for the Laplace equations which can be solved 
explicitly. The fact that the complex Ernst equation which takes the role of the Laplace 
equation in the relativistic case is completely integrable gives rise to the hope that bound- 
ary value problems might be solvable here at least in selected cases. The unifying frame- 
work for both the Laplace and the Ernst equation is provided by methods from soliton 
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theory, so-called Riemann-Hilbert problems: the scalar problem for the Laplace equation 
can be always solved with the help of a generalization of the Cauchy integral (see [16j and 
references given therein), a procedure which leads to the Poisson integral for distributional 
densities. Choosing the contour of the Riemann-Hilbert problem appropriately one can 
construct solutions to the Laplace equation which are everywhere regular except at a disk 
where the function is not differentiable. Similarly one can treat the relativistic case where 
the matrix Riemann-Hilbert problem can be related to a linear integral equation. It was 
shown in |T7| that the matrix problem for the Ernst equation can be always gauge trans- 
formed to a scalar problem on a Riemann surface which can be solved explicitly in terms 



of Korotkin's finite gap solutions |L8j for rational Riemann-Hilbert data. In this sense 
these solutions can be viewed as a generalization of the Poisson integral to the relativistic 
case. 

Whereas the Poisson integral contains one free function which is sufficient to solve bound- 
ary value problems for the scalar gravitational potential, the finite gap solutions contain 
one free function and a set of complex parameters, the branch points of the Riemann 
surface. Thus one cannot hope to solve general boundary value problems for the complex 
Ernst potential within this class because this would imply the choice to specify two free 
functions in the solution according to the boundary data. This means that one can only 
solve certain classes of boundary value problems on a given compact Riemann surface. 
In the first article we investigate the implications of the underlying Riemann surface on 
the multipole moments and the boundary values taken at a given boundary. The rela- 
tions will be given for general genus of the surface and will be discussed in detail in the 
case of genus 1 (elliptic surface) and genus 2, which is the simplest case with generic 
equatorial symmetry. It is shown that the solution of boundary value problems leads in 
general to non-linear integral equations. We can identify however classes of boundary 
data where only one linear integral equation has to be solved. Special attention will be 
paid to counter-rotating dust disks which will lead us to the construction of the solution 
for constant angular velocity and constant relative density which was presented in ]Tj|. 
It contains as limiting cases the static solutions of Morgan and Morgan [B(J and the disk 



with only one matter stream by Neugebauer and Meinel ||21|| . The potentials of the re- 
sulting spacetime at the axis and the disk are presented in the second article, the physical 
features as the ultrarelativistic limit, the formation of ergospheres, multipole moments 
and the energy-momentum tensor are discussed in the third article. 
The present article is organized as follows. In section ^| we discuss Newtonian dust disks 
with Riemann-Hilbert methods and relate the corresponding boundary value problems 
to an Abelian integral equation. The relativistic field equations and the boundary con- 
ditions for counter-rotating dust disks are summarized in section || Important facts on 
hyperelliptic Riemann surfaces which will be used to discuss Korotkin's class of solutions 
to the Ernst equation are collected in section |I[ In section [5], we establish relations for 
the corresponding Ernst potentials on the axis on a given Riemann- surface of arbitrary 
genus. The found relation limits the possible choice of the multipole moments. We dis- 
cuss in detail the elliptic and the genus 2 case with equatorial symmetry. This analysis 
is extended to the whole spacetime in section || which leads to a set of differential and 
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algebraic equations which is again discussed in detail for genus 1 and 2. The equations 
for genus 2 are used to study differentially counter-rotating dust disks in section |7|: We 
discuss the Newtonian limit of disks of genus 2. As a first application of this constructive 
approach we derive the class of counter-rotating dust disks with constant angular velocity 



and constant relative density of [19|. We prove the regularity of the solution up to the 



ultrarelativistic limit in the whole spacetime except the disk and conclude in section 



2 Newtonian dust disks 

To illustrate the basic concepts used in the following sections, we will briefly recall some 
facts on Newtonian dust disks. In Newtonian theory, gravitation is described by a scalar 
potential U which is a solution to the Laplace equation in the vacuum region. We use 
cylindrical coordinates p, ( and <p and place the disk made up of a pressureless two- 
dimensional ideal fluid with radius po in the equatorial plane ( = 0. In Newtonian theory 
stationary perfect fluid solutions and thus also the here considered disks are known to be 
equatorially symmetric. 

Since we concentrate on dust disks, i.e. pressureless matter, the only force to compensate 
gravitational attraction in the disk is the centrifugal force. This leads in the disk to (here 
and in the following f x = |^ ) 

u p = n 2 (p) P , (i) 

where fl(p) is the angular velocity of the dust at radius p. Since all terms in ([!]) are 
quadratic in Q there are no effects due to the sign of the angular velocity. The absence 
of these so-called gravitomagnetic effects in Newtonian theory implies that disks with 
counter-rotating components will behave with respect to gravity exactly as disks which 
are made up of only one component. We will therefore only consider the case of one 
component in this section. Integrating ([!]) we get the boundary data U (p, 0) with an 
integration constant Uq = U(0, 0) which is related to the central redshift in the relativistic 
case. 

To find the Newtonian solution for a given rotation law Q{p), we thus have to construct a 
solution to the Laplace equation which is everywhere regular except at the disk where it 
has to take on the boundary data (|1]). At the disk the normal derivatives of the potential 
will have a jump since the disk is a surface layer. Notice that one only has to solve the 
vacuum equations since the two-dimensional matter distribution merely leads to boundary 
conditions for the Laplace equation. In the Newtonian setting one thus has to determine 
the density for a given rotation law or vice versa, a well known problem (see e.g. and 
references therein) for Newtonian dust disks. 

The method we outline here has the advantage that it can be generalized to some extent 
to the relativistic case. We put po = 1 without loss of generality (we are only considering 
disks of finite non-zero radius) and obtain U as the solution of a Riemann-Hilbert problem 



(see e.g. |I6[ and references given therein), 
Theorem 2.1: 

Let InG 6 C 1,a (T) and T be the covering of the imaginary axis in the upper sheet of S c 
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between — i and i where Eo is the Riemann surface of genus given by the algebraic relation 
/Zq(t) = (r — () 2 + p 2 . The function G has to be subject to the conditions G(f) = G{r) 
and G(—t) = G(r). Then 

U( P ,C) = -—[ ' nG{T)dT (2) 
4th Jr . / (r _ C) 2 + p 2 y } 



is a real, equatorially symmetric solution to the Laplace equation which is everywhere 
regular except at the disk ( = 0, p < 1. The function InG is determined by the boundary 
data U(p, 0) or the energy density a of the dust (2ira = in units where the velocity of 
light and the Newtonian gravitational constant are equal to 1) via 

or 

lnGW = 4 / 1 7FM' (4) 

respectively where t = —it. 

The occurrence of the logarithm in (§) is due to the Riemann-Hilbert problem with the 
help of which the solution to the Laplace equation was constructed. We briefly outline 
the 

Proof: 

It may be checked by direct calculation that U in (H) is a solution to the Laplace equation 
except at the disk. The reality condition on G leads to a real potential, whereas the 
symmetry condition with respect to the involution r — >• — r leads to equatorial symmetry. 
At the disk the potential takes due to the equatorial symmetry the boundary values 



1 fp lnG(t) 

and 



1 fi dtQnGW 

Both equations constitute integral equations for the 'jump data' InG of the Riemann- 
Hilbert problem if the respective left-hand side is known. The equations @ and (jg) are 
both Abelian integral equations and can be solved in terms of quadratures, i.e. (||) and 
(|). To show the regularity of the potential U we prove that the integral (§) is identical 
to the Poisson integral for a distributional density which reads at the disk 

u[p) = -2 /' r j* = -4 /' wMmtm, ( 7 ) 

yJ(p + P ') 2 -4pp' cos (f) J o P + P 

where k(p, p') = 2^/p/7 /(p + p') and where K is the complete elliptic integral of the first 
kind. Eliminating In G in (H) via (f|) we obtain after interchange of the order of integration 



U = --\ U c H -K[ H -\dp'+ U C K ^)dp'\ (8) 
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which is identical to (0) since K(2Vk/(l+k)) = (l+k)K(k). Thus the integral (g) has the 
properties known from the Poisson integral: it is a solution to the Laplace equation which 
is everywhere regular except at the disk where the normal derivatives are discontinuous. 
This completes the proof. 

Remark: We note that it is possible in the Newtonian case to solve the boundary value 
problem purely locally at the disk. The regularity properties of the Poisson integral then 
ensure global regularity of the solution except at the disk. Such a purely local treatment 
will not be possible in the relativistic case. 

The above considerations make it clear that one cannot prescribe both U at the disk (and 
thus the rotation law) and the density independently. This just reflects the fact that the 
Laplace equation is an elliptic equation for which Cauchy problems are ill-posed. If In G is 
determined by either (§) or @ for given rotation law or density, expression (0) gives the 
analytic continuation of the boundary data to the whole spacetime. In case we prescribe 
the angular velocity, the constant Uq is determined by the condition In (7(1) = which 
excludes a ring singularity at the rim of the disk. For rigid rotation (Q = const), we get 
e.g. 

lnG(r) = 4ft 2 (r 2 + l) (9) 
which leads with (|2[) to the well-known Maclaurin disk. 



3 Relativistic equations and boundary conditions 

It is well known (see |[22| ) that the metric of stationary axisymmetric vacuum spacetimes 
can be written in the Weyl-Lewis-Papapetrou form 

ds 2 = -e 2U (dt + ad0) 2 + e~ 2U (e 2k (dp 2 + d( 2 ) + p 2 d0 2 ) (10) 

where p and ( are Weyl's canonical coordinates and d t and are the two commuting 
asymptotically timelike respectively spacelike Killing vectors. 

In this case the vacuum field equations are equivalent to the Ernst equation for the 
complex potential / where / = e 2U + ib, and where the real function b is related to the 
metric functions via 

b z = --e 4 V- (11) 
P 

Here the complex variable z stands for z — p + iC With these settings, the Ernst equation 
reads 

1 2 

fzz ~\~ 777 ; Z\\fz ~l" fz) ~~% : jfzfz j (12) 
2(z + z) f + f 

where a bar denotes complex conjugation in C. With a solution /, the metric function U 
follows directly from the definition of the Ernst potential whereas a can be obtained from 
( p|) via quadratures. The metric function k can be calculated from the relation 

k z = 2p(U z ) 2 -^e 4U (a z ) 2 . (13) 
2p 
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The integrability condition of ([TTJ) and ( |T3"D is the Ernst equation. For real /, the Ernst 
equation reduces to the Laplace equation for the potential U. The corresponding solutions 
are static and belong to the Weyl class. Hence static disks like the counter-rotating disks 
of Morgan and Morgan [^(J can be treated in the same way as the Newtonian disks in the 



previous section. 

Since the Ernst equation is an elliptic partial differential equation, one has to pose bound- 
ary value problems. The boundary data arise from a solution of the Einstein equations 
in the matter region. In our case this will be an infinitesimally thin disk made up of two 
components of pressureless matter which are counter-rotating. These models are simple 
enough that explicit solutions can be constructed, and they show typical features of gen- 
eral boundary value problems one might consider in the context of the Ernst equation. 
It is also possible to study explicitly the transition from a stationary to a static space- 
time with a matter source of finite extension for these models. Counter-rotating disks of 



infinite extension but finite mass were treated in and p3 |, disks producing the Kerr 



metric and other metrics in flTI 



To obtain the boundary conditions at a relativistic dust disk, it seems best to use Israel's 
invariant junction conditions for matching spacetimes across non-null hypersurfaces . 
Again we place the disk in the equatorial plane and match the regions V ± (±£ > 0) 



at the equatorial plane. This is possible with the coordinates of (|10| ) since we are only 
considering dust i.e. vanishing radial stresses in the disk. The jump 7 Q/ g = — K~p in 
the extrinsic curvature K a p of the hypersurface ( = with respect to its embeddings into 
V ± = {±( > 0} is due to the energy momentum tensor S a p of the disk via 

- 8nS af3 = ^ - h a ^\ (14) 

where h is the metric on the hypersurface (greek indices take the values 0, 1, 3 corre- 
sponding to the coordinates t, p, 0). As a consequence of the field equations the energy 
momentum tensor is divergence free, S°f = where the semicolon denotes the covariant 
derivative with respect to h. 

The energy-momentum tensor of the disk is written in the form 

= a+u^ul + a^u1u u _, (15) 

where the vectors linear combination of the Killing vectors, {u±) = (1, 0, ±f2(p)). 

This has to be considered as an algebraic definition of the tensor components. Since 
the vectors u± are not normalized, the quantities <j± have no direct physical significance, 
they are just used to parametrize S^ u . The energy-momentum tensor was chosen in 
a way to interpolate continuously between the static case and the one-component case 
with constant angular velocity. An energy-momentum tensor with three independent 
components can always be written as 

= a* p v"v v + p* p w "vf, (16) 

where v and w are the unit timelike respectively spacelike vectors (t> M ) = iVi(l, 0, uo^) 
and where (u> M ) = N 2 (k, 0, 1). This corresponds to the introduction of observers (called 
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isotropic observers (FIOs) in [IT]) for which the energy- momentum tensor is diagonal. 



The condition w^v M = determines n in terms of tu^ and the metric, 

K = — . (17) 

<?oo + w<f>go3 

If Pp/<Jp < 1 the matter in the disk can be interpreted as in jZIJ either as having a purely 
azimuthal pressure or as being made up of two counter-rotating streams of pressureless 
matter with proper surface energy density cx*/2 which are counter-rotating with the same 



angular velocity yPp/<J* } 

= \*{U'lU u + + mui) (18) 

where (U±) = U*(v^ ± \jp* v j 'o*w^) is a unit timelike vector. We will always adopt the 
latter interpretation if the condition p*/cr* < 1 is satisfied which is the case in the example 
we will discuss in more detail in section 7. The energy-momentum tensor (|TJ|) is just the 
sum of two energy-momentum tensors for dust. Furthermore it can be shown that the 
vectors U± are geodesic vectors with respect to the inner geometry of the disk: this 
is a consequence of the equation S^f = together with the fact that U± is a linear 
combination of the Killing vectors. In the discussion of the physical properties of the disk 
we will refer only to the measurable quantities lu^, a* and p* which are obtained by the 
introduction of the FIOs whereas a± and Q are just used to generate a sufficiently general 
energy-momentum tensor. 

To establish the boundary conditions implied by the energy-momentum tensor, we use 
Israel's formalism f24|. Equation S°f = leads to the condition 



U p (l + 2 7 fia + ftV) + fia p ( 7 + fta) + n 2 p(pU p - l)e~ iU = 0, (19) 

where 

7(P) = — pr- (20) 

The function j(p) is a measure for the relative energy density of the counter-rotating 
matter streams. For 7 = 1, there is only one component of matter, for 7 = 0, the matter 
streams have identical density which leads to a static spacetime of the Morgan and Morgan 
class. 

As in the Newtonian case, one cannot prescribe both the proper energy densities <j± and 
the rotation law Q at the disk since the Ernst equation is an elliptic equation. For the 
matter model (|15|), we get at the disk 
Theorem 3.1: 

Let a(p) = cr + (p) + cr_(p) and let R(p) and S(p) be given by 

R=(a + l)e^, (2!) 

and 

*M-^- « 22 > 
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Then for prescribed Q(p) and S(p), the boundary data at the disk take the form 

k ~~ l 2Rp fp+ R 6 ■ (23) 



Let a be given by a = ere . Then for given density a and 7, the boundary data read, 

(P 2 + ((e 2U ) p (e 2U ) ( + b p b c ) 2 - 2 P e 2U {e 2U ) ( ((e 2U ) p (e^) f + b p b c ) + b 2 / u = 0, 

(24) 

and 



b p - a 



( e2U ) P ( e2U ) ( + & A)) 2 + STrpae'W ({e 2U ) p (e 2C7 ) ( + b p b c ) = 0. (25) 



Proof: 

The relations fli"4}) lead to 



- 47re (fc - l/) S'oo = (Jfe C - 2C/ C ) e 2C/ 
-47re( fe - c/ )( ( So3-a^ Q o) = -Ja c e 2C/ , 



2 

-4^-^(533- 20^3 + a 2 S , 00 ) = -k c p 2 e- 2U , (26) 



where 



S 00 = ae w (l + Q 2 a 2 + 2fia 7 ) , 

S03 - O'S'oo = -vp 2 Q (Qa + 7) , 
S 33 - 2aS , 03 + a 2 ^00 = c^Ve" 4 ". (27) 



One can substitute one of the above equations by ( |19| ) in the same way as one replaces 
one of the field equations by the covariant conservation of the energy momentum tensor 
in the case of three-dimensional ideal fluids. This makes it possible to eliminate from 
(p6|) and to treat the boundary value problem purely on the level of the Ernst equation. 
The function k will then be determined via flHf ) with the found solution of the Ernst 
equation. It is straight forward to check the consistency of this approach with the help of 

(0)- 

If Q and 7 (and thus S) are given, one has to eliminate a from ( p6|) and (p7). This can 
be combined with flTJ) and (|TTp to (|23]). 

If the function 7 and a are prescribed (this makes it possible to treat the problem com- 
pletely on the level of the Ernst equation), one has to eliminate Q from (|19D, ( p6"D and 
( p?D which leads to (|24T ) and (p5|). This completes the proof. 

Remark: For given fi(p) and S(p), equation (|i~9"D is an ordinary non-linear differential 
equation for e 2U , 

(R 2 - p 2 ) p e 2U - 2Re AU (1) = (R 2 - p 2 - 5e w ) (e 2U ) . (28) 
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For constant f2 and 7 we get 

i? 2_ p 2 + 5e 4C/ = 2 2 „ (2g) 

A 

where A = 2tt 2 e- 2Uo . 

For given boundary values as in Theorem 3.1, the task is to to find a solution to the Ernst 
equation which is regular in the whole spacetime except at the disk where it has to satisfy 
two real boundary conditions. In the following we will concentrate on the case where the 
angular velocity Q and the relative density 7 are prescribed. 



4 Solutions on hyperelliptic Riemann surfaces 

The remarkable feature of the Ernst equation is that it is completely integrable which 
means that the Riemann- Hilbert techniques used in the Newtonian case can be applied 
here, too. This time, however, one has to solve a matrix problem (see e.g. JlTj and 



references given therein) which cannot be solved generally in closed form. In |T7| it was 
shown that the problem can be gauge transformed to a scalar problem on a four-sheeted 
Riemann surface. In the case of rational 'jump data' of the Riemann-Hilbert problem, this 
surface is compact and the corresponding solutions to the Ernst equation are Korotkin's 
finite gap solutions [13]. In the following we will concentrate on this class of solutions and 



investigate its properties with respect to the solution of boundary value problems. 

4.1 Theta functions on hyperelliptic Riemann surfaces 

We will first summarize some basic facts on hyperelliptic Riemann surfaces which we will 
need in the following. We consider surfaces S of genus g which are given by the rela- 
tion n 2 (K) = (K + iz)(K - iz) U 9 i=1 (K - Ei)(K - Ei) where the Ei do not depend on 
the physical coordinates z and z. We introduce the standard quantities associated with 
a Riemann surface (see [p5H ), with respect to the cut system of figure 1 (we order the 
branch points with ImEi < in a way that ReEi < Re-E^ < . . . < KeE g and assume for 
simplicity that the real parts of the Ei are all different; we write Ei = on + the g nor- 
malized differentials of the first kind du>i defined by § a . dujj = 2ni5ij, and with Pq = —iz 
the Abel map u>i(P) = Jp du>i which is defined uniquely up to periods. Furthermore, 
we define the Riemann matrix II with the elements = § b . duj, and the theta func- 

tion e[m](z) = }- iVeZgeX p{l(n(iV+^),(iV + ^)) + ((^ + 7rim 2 ),(iV+f))} with 

half integer characteristic [m] = ™ 2 and mj,m 2 = 0,1 ((N,z) = Yh=i NiZi) . A charac- 
teristic is called odd if (m 1 ,^. 2 ) 7^ mod 2. The normalized (all a-periods zero) differ- 
ential of the third kind with poles at Pi and Pi and residue +1 and —1 respectively will 
be denoted by dup 1 p 2 . A point P e E will be denoted by P = (K, ±fj,(K)) or K ± (the 
sheets will be defined in the vicinity of a given point on S, e.g. 00). 
The theta functions are subject to a number of addition theorems. We will need the 
ternary addition theorem which can be cast in the form 
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Figure 1: Homology basis of E 



Theorem 4.1: Ternary addition theorem 

Let [mi] = [m},mf] (i — 1, . . . ,4) be arbitrary real 2g- dimensional vectors. Then 



9[mi](u + v)Q[m 2 ]{u - v)e[m 3 ]{0)Q[m 4 }(0) (30) 
= — exp(— 47ri(m\, a 2 ))Q[ni + a](u)Q[ri2 + a](u)Q[n 3 + a](i7)9[n 4 + a](v), 

29 2ae(Z 2 ) 2 f 

where a = (a 1 , a 2 ), and (m 1; . . . , 777.4) = (774, . . . , n4)T wit/i 
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T=- n : - : . pi) 



n k P k with Pi e S 



Each 1 in T denotes the g x g identity matrix. 
For a proof see e.g. P5[ . 

Let us recall that a divisor X on E is a formal symbol X = 774P1 + . 
and rij 6 Z. The degree of a divisor is X^=i n i- The Riemann vector i^^ is defined by the 
condition that 0(u(W) +K R ) = if W is a divisor of degree g — 1 or less. We use here and 
in the following the notation u(W) = fp^du = Y%=i w{wi). We note that the Riemann 
vector can be expressed through half-periods in the case of a hyperelliptic surface. 
The quotient of two theta functions with the same argument but different characteristic 
is a so-called root function which means that its square is a function on S. One can prove 
(see [^6| and references therein) 
Theorem 4.2: Root functions 

Let Qi, % = 1, . . . , 2g + 2, be the branch points of a hyperelliptic Riemann surface S g of 
genus g and Aj = u>(Qj) with u>(Qi) = 0. Furthermore let {ii, . . . ,i g } and {ji, . . . ,j g } 
be two sets of numbers in {1, 2, . . . , 2g + 2}. Then the following equality holds for an 
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arbitrary point P 6 S 9 , 



e[K R + n=iA k ]^(P)) _ (K-E ll )...(K-E lg ) 
6 [K R + ELi A jh ] (u(P)) Cl \ (K -E h )...{K- E jg \ 



•where c\ is a constant independent on K . Let X = P\ + . . . + P g with Pj 
a divisor of degree g on E 9 then the following identity exists, 



(32) 



Q[K R + Ai] (cu(X)) 



Q[K R + Aj](u(X)) ~ 2 t=i\ {Kk-Q. 



(33) 



where C2 is a constant independent on the Kk. 

The notion of divisors makes it possible to state Jacobi's inversion theorem in a very 
compact form, 

Theorem 4.3: Jacobi inversion theorem 

Let A, B G S be divisors of degree g and u G C g . Then for given B and u, the equation 
uj(A) — u(B) = u for the divisor A is always solvable. 

For a proof we refer the reader to the standard literature, e.g. P5 |. We remark that the 
divisor may not be uniquely defined in the general case which means that one or more 
Pi G A can be freely chosen. We will not consider such special cases in the following and 



refer the reader for the so-called special divisors to the literature as |p6 
For divisors A — B of degree zero, one can formulate Abel's theorem. 
Theorem 4.4: Abel's theorem 

Let A,B G E be divisors of degree n subject to the relation u(A) — u)(B) = 0. Then A 
and B are the set of zeros respectively poles of a meromorphic function F . 
For a proof see [p5|| . We remark that this function is a rational function on the surface 
cut along the homology basis. We have the 
Corollary 4.5: 

Let the condition of Abel's theorem hold. Then the following identity holds for the integral 
of the third kind 

F(P) 



B 



PQ 



In 



F(QY 



(34) 



4.2 Solutions to the Ernst equation 

We are now able to write down a class of solutions to the Ernst equation on the surface 
S. 

Theorem 4.5: 

Let the Riemann surface S be given by the relation fi 2 (K) = (K + iz)(K — \z) Y[ g i=l (K — 
Ei)(K — Ei), let u be the vector with the components Ui — J T InGdui where T is as in 
theorem 2.1, let G be subject to the condition G(t) = G(f), and let [m] = [m l ,m 2 ] with 
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m\ = and m 2 arbitrary for % = 1, . . . ,g be a theta characteristic. Then the function f 
given by 

ft r\ ®W(u(oo + )+u) \ 1 f \ 
f{P, = nr T7 7 _x , — r exp I — / In G r duoo+oo- r \ , 35 
6|mJ(u;(oo )+«) [27riJ J 

zs a solution to the Ernst equation. 

This class of solutions was first given by Korotkin [18|, the straight forward continuous 
limit leading to the above form can be found in [P7|, |28fl . For the relation to Riemann- 
Hilbert problems see |L7|]. In the case genus 0, the Ernst potential is real, and we get a 



solution of the Weyl class in the form (g). For higher genus, these solutions are in general 
non-static and thus generalize (0) to the stationary case. 



In [£9| £(]] it was possible to identify a physically interesting subclass. 
Theorem 4.6: 

Let the conditions of Theorem 4-5 hold, and in addition let £ be a hyperelliptic Riemann 
surface of even genus g = 2n given by fi 2 (K) — (K + iz)(K — iz) TYt =1 (K 2 — Ef)(K 2 — 
Ef), let the function G be subject to the condition G(—r) = G(t), and let [n] be the 
characteristic with n\ = and nf = 1 . Then the function f given by 

e[n](uj(oo + ) +u) / 1 , . 



r 



is an equatorially symmetric solution to the Ernst equation (f(—C) — f(C)) which is 
everywhere regular except at the disk if ©(u;(oo~) + u) ^ 0. 

For a proof see |29], where one can also find how the characteristic can be generalized. 



In the following we will only use the characteristic of the above theorem. 

A quantity of special interest is the metric function a. In |18] it was shown that one can 



relate it directly to theta functions without having to perform the integration of (|TT|), 



2u - , 6(^(oo ) +u) 



Z:=(a- a )e 2U = D M - In \ \ ' ' (37) 

(3[n\(uj(oo ) + u) 

where DpF(u(P)) denotes the coefficient of the linear term in the expansion of the func- 
tion F(u(P)) in the local parameter in the neighborhood of P, where is the Riemann 
theta function with the characteristic [m] and m\ = mf = 0, and where the constant ao 
is defined by the condition that a vanishes on the regular part of the axis. 
It is possible to give an algebraic representation of the solutions (|36|) (see |3lJ and [p2| ). 
We define the divisor X = J2f=i Ki as the solution of the Jacobi inversion problem (i = 

Ui(X) - Ui(D) = — / lnG^^ =: u h (38) 
where the divisor D = J2f=i Ei- With the help of these divisors, we can write (|36| ) in the 



form 

, „ f x tHt 1 r , tHt , . 
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Since the U{ in are just the periods of the second integral in Q3TS|), they are subject to a 
system of differential equations, the so-called Picard-Fuchs system (see |3C| and references 
given therein). In our case this leads to 

^ (^-Po)^ Q; - = 2 (40) 

and 

(ln/), = E ( ^ w "^ 0) 5 r1 ^. (41) 
Solving for the i^ n , z , n = 1, . . . , g, we get 

^ 2 = L (42) 

Additional information follows from the reality of the Uj which implies w(X) — u)(D) = 
u>(X) — uj{D). Using Abel's theorem on this condition, we obtain the relation for an 
arbitrary K G C 



1 - x 2 ) f[(K - Ki)(K - Kt) = f[(K - Ei)(K - E,) - (K - P )(K - P )Q 2 2 (K), (43) 



f[(K-K l )(K-K t ) = f[ 

%=i i=i 

where with purely imaginary Xi, x 

Q 2 {K) =x + x x K + ... + x g . 2 K 9 ~ 2 + xK 9 ~\ (44) 



Since (^3j) has to hold for all K G C, it is equivalent to 2g real algebraic equations for the 
Ki if the Xi are given. Using ( |34]) and ( P5| ) we find 

4 - ^ (45) 

which implies x = ibe~ 2U . 

Remark: To solve boundary value problems with the class of solutions (|36D , one has 
two kinds of freedom: the function G as before and the branch points E^ of the Riemann 
surface as a discrete degree of freedom. Since one would need to specify two free functions 
to solve a general boundary value problem for the Ernst equation, it is obvious that one 
can only solve a restricted class of problems on a given surface, and that one cannot expect 
to solve general problems on a surface of finite genus. But once one has constructed a 
solution which takes the imposed boundary data at the disk, one has to check the condition 
Q(u(oo~)+u) ^ in the whole spacetime to actually prove that one has found the desired 
solution: a solution that is everywhere regular except at the disk where it has to take the 
imposed boundary conditions. 

There are in principle two ways of generalizing the approach used for the Newtonian case: 
One can eliminate O from the two real equations (|23|) and enter the resulting equation 
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with a solution ( ]56D on a chosen Riemann surface. This will lead for given 7 to a non- 
linear integral equation for In G. In general there is little hope to get explicit solutions to 
this equation (for a numerical treatment of differentially rotating disks along this line in 
the genus 2 case see [Q). Once a function G is found, one can read off the rotation law 
Q on a given Riemann surface from fl36l) . Another approach is to establish the relations 
between the real and the imaginary part of the Ernst potential which exist on a given 
Riemann surface for arbitrary G. The simplest example for such a relation is provided by 
the function w = e 1 ^ which is a function on a Riemann surface of genus 0, where we have 
obviously \w\ = 1. As we will point out in the following, similar relations also exist for an 
Ernst potential of the form fl36l), but they will lead to a system of differential equations. 
Once one has established these relations for a given Riemann surface, one can determine 
in principle which boundary value problems can be solved there (in our example which 
classes of functions Q, 7 can occur) by the condition that one of the boundary conditions 
must be identically satisfied. The second equation will then be used to determine G as 
the solution of an integral equation which is possibly non-linear. Following the second 
approach, we want to study the implications of the hyperelliptic Riemann surface for the 
physical properties of the solutions. 



5 Axis Relations 

In order to establish relations between the real and the imaginary part of the Ernst 
potential, we will first consider the axis of symmetry (p = 0) where the situation simplifies 
decisively. In addition the axis is of interest since the asymptotically defined multipole 
moments fl34| , |35] can be read off there |36] . 



On the axis the Ernst potential can be expressed through functions defined on the Rie- 
mann surface £' given by // 2 = nf = i(K — Ei)(K — Ei), i.e. the Riemann surface obtained 
from E by removing the cut [Po, Pq] which just collapses on the axis. We use the notation 
of the previous section and let a prime denote that the corresponding quantity is defined 
on the surface The cut system is as in the previous section with [Ex, Ex] taking the role 
of [P ,-Po] ( a h fr-cuts cross [Ex, Ex])- We choose the Abel map in a way that u'(Ex) = 0. 
It was shown in [3(J that for genus g > 1 the Ernst potential takes the form (for ( > 0) 



(m do/ + u') - exp - K 00+ + u g ))0 (f£ du' + u') T ^ 

mO = -p- + { 3 P . fe 7+? \ (46) 

■& (/ c °f dcu' - u>) - exp(-(^(oo+) - u g ))$ (/ c °° dco' - u>) 

where d is the theta function on £' with the characteristic a' { = 0, f3[ = 1 for i = 1, . . . , g — 

1, where / = ±? f In G(r)du' 00+00 - (r), du g = du; c - c +, and where u g = =jjr JlnG(r)da; g (r). 

r r 
Notice that the u[ and / are constant with respect to (. The only term dependent both on 

G and on ( is u g . To establish a relation on the axis between the real and the imaginary 

part of the Ernst potential independent of G, the first step must be thus to eliminate u g . 

We can state the following 
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Theorem 4.1: 

The Ernst potential faffl ) satisfies for g > 1 the relation 

Pi(0ff + P2(()b + P 3 (0 = 0, (47) 

where the Pi are real polynomials in ( with coefficients depending on the branch points 
Ei and the g real constants J r In Gr l dr / //(r) with i = 0, . . . , g — 1. The degree of the 
polynomials P\ and P3 is 2g — 3 or less, the degree of P2 is 2g — 2 or less. 
To prove this theorem we need the fact that one can express integrals of the third kind 
via theta functions with odd characteristic denoted by i? , 

exp(-. 9(o o+)) = - ^U-^—^Jj . (48) 

Proof: 

The first step is to establish the relation 

Aff + Bib + 1 = 0, (49) 

where 

and 
Be T = 

1? [u' + / c °f ' du/J 1? (u' + J c + ' cL/J 

(51) 

which may be checked with (^) by direct calculation. The reality properties of the 
Riemann surface £' and the function G imply that A is real and that B is purely imaginary. 
We use the addition theorem fl3~0|) with [mi] = . . . = [mj equal to the characteristic of 1? 
for fl5Ti| ) to get 

£ exp(-47ri(m},a 2 ))^ 2 [a]( M , + ^ / (oo-))^ 2 [a](^ / (C + )) 

£ exp(-47ri(m{,a 2 ))^ 2 [a](M' + w / (oo+))^ 2 [a](cu / (C+))' 1 J 

2ae(Z 2 ) 2 9 

This term is already in the desired form. Using the relation for root functions fl32|), one 
can directly see that the right-hand side is a quotient of polynomials of order g — 1 or 
lower in £. For (|5"1"D we use (|48| ) with [mj = [777,2] = [i£#] as the characteristic of the odd 
theta function d and let [m 3 ] = [7714] be equal to the characteristic of d. The addition 
theorem (I3TI) then leads to 



Ae2 , = tfftf + /~-da/)tf (tf + ff- do/) 
(u' + /~ + do)') 1? (u' + /+ dw') 



e -^(oo+)^ (y + joo+ ^ 


ll?( 


u' + do;') 




[vf + J c °f do/) 




(1/' + /«" do;') 


1?l 


[u> + ; c °? + <l/ 


) * («' + 


do/) 


1 



E exp(-47ri(mi,a 2 ))e /2 [n + a](M / + w / (oo-))7? 2 [a](cu / (C + )) 

d J 2ag(Z 2 )2g 

6 E exp(-47ri(ml,a 2 ))i? 2 [a](M / + ^ / (oo+))i? 2 [a](u; / (C + )) 

2aG(Z 2 ) 2 9 
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2_, exp(-47ri(m^ a 2 )) i!^ x 



2ae(Z 2 ) 2 9 



i? 2 (0) 



^ 2 [a](,/(oo+)-a/(C + )) , ^[a](a/(cx)+)-a/(C-)) \ , . 

, ^'(00+) - u/(C+)) + ^(w'(oo+)-a;'(C-)) J' 1 J 



where n follows from rh as in Theorem 4.1. The first fraction in ([53]) is again the quotient 
of polynomials of degree g — 1 in £ for the same reasons as above. But since the quotient 
must vanish for £ — > oo, the leading terms in the numerator just cancel. It is thus a 
quotient of polynomials of degree g — 2 or less in the numerator and g — 1 or less in the 
denominator. To deal with the quotients $ 2 [a](u/(oo + ) — c</(£ ± ))/'# 2 (u/(oo + ) — dj / (^ ± )), we 
define the divisors = Tf 1 + . . . + TgL 1 as the solutions of the Jacobi inversion problems 
0/(1*) - a/(r) = lu'(oo + ) + w'(C ± ) where K is the divisor Y = E 1 + ... + E g _ x . Abel's 
theorem then implies for arbitrary K E C 

fiV - Tt){K -Q = (K- A±f nV - Et) -(K- E g ) f[(K- E,), (54) 

i=l i=l i=l 

where 

±= — /AO — (55) 

Let Qj be given by the condition [Qj + If^] = [a], i.e. Q., is a branch point of Then 
we get for the quotient 

m^t^gn = cmst g gz* (56) 



where const is a ^-independent constant. With the help of (|54|), it is straight forward to 
see that for Qj 6 F, the theta quotient is just proportional to (£ — Ex)/(( — Qj) whereas 
for ^ y, the term is proportional to (£ — E 1 )(Qj — A ± ) 2 /(C — Qj). Using ( |55|) one 
recognizes that the terms containing roots just cancel in (|53~D . The remaining terms are 
just quotients of polynomials in ( with maximal degree g in the numerator and g — 2 in 
denominator. This completes the proof. 

Remark: The remaining dependence on G through v! and I can only be eliminated 
by differentiating relation ([46]) g times. If we prescribe e.g. the function b on a given 
Riemann surface (this just reflects the fact that the function G can be freely chosen in 
( pj|) ), we can read off e 2U from (|4"7]). To fix the constants related to G in ([4"7]) one needs 
to know the Ernst potential and g — 1 derivatives at some point on the axis where the 
Ernst potential is regular, e.g. at the origin or at infinity, where one has to prescribe the 
multipole moments. If the Ernst potential were known on some regular part of the axis, 
one could use ([IT]) to read off the Riemann surface (genus and branch points). Equation 
([46|) is then an integral equation for G for known sources. This just reflects a result of 
|37J that the Ernst potential for known sources can be constructed via Riemann- Hilbert 



techniques if it is known on some regular part of the axis. 
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In practice it is difficult to express the coefficients in the polynomials Pi via the constants 
u[ and /, and it will be difficult to get explicit expressions. We will therefore concentrate 
on the general structure of the relation fl47f) , its implications on the multipoles and some 
instructive examples. Let us first consider the case genus 1 which is not generically 
equator ially symmetric. In this case the Riemann surface £' is of genus 0. One can use 
formula for the axis potential if one replaces the theta functions by 1. We thus end 
up with 

ff-2b ( —^e I = e 21 . (57) 

Pi 

Here the only remaining G-dependence is in /. If f = /(0,0) is given, e 1 follows from 
f f +2bodie I I (3\ = e 21 , if the in general non-real mass M is known, the constant e 1 follows 
from 1 + 4ImMe 7 //3i = e 21 . In the latter case the imaginary part of the Arnowitt-Deser- 
Misner mass (this corresponds to a NUT-parameter) will be sufficient. Differentiating 



(p7|) once will lead to a differential relation between the real and the imaginary part of 
the Ernst potential which holds for all G, which means it reflects only the impact of the 
underlying Riemann surface on the structure of the solution. 

Remark: For equatorially symmetric solutions, one has on the positive axis the relation 
/(— C)/(C) = 1 ( see HE P9| ) this is to be understood in the following way: the function |£| 
is even in (, but restricted to positive ( it seems to be an odd function, and it is exactly 
this behavior which is addressed by the above formula). This leads to the conditions 

P 1 (-() = -P 3 ((), P 2 (-0 = P 2 (C). (58) 

The coefficients in the polynomials depend on the g/2 integrals J r dr In Gt 2i /h'(t) (i = 

0. . . . , g/2 — 1 and the branch points. 

The simplest interesting example is genus 2, where we get with E\ = a + i/3 

ff(C - Ci) + ^(C 2 -a- C 2 2 )b = C + Ci, (59) 

1. e. a relation which contains two real constants C±, C2 related to G. In case that the 
Ernst potential at the origin is known, one can express these constants via /o- A relation 
of this type, which is as shown typical for the whole class of solutions, was observed in 



the first paper of [21] for the rigidly rotating dust disk. 



6 Differential relations in the whole spacetime 

The considerations on the axis have shown that it is possible there to obtain relations 
between the real and the imaginary part of the Ernst potential which are independent 
of the function G and thus reflect only properties of the underlying Riemann surface. 
The found algebraic relations contain however g real constants related to the function 
G, which means that one has to differentiate g times to get a differential relation which 
is completely free of the function G. These constants were just the integrals u' and / 
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which are only constant with respect to the physical coordinates on the axis where the 
Riemann surface S degenerates. Thus one cannot hope to get an algebraic relation in the 
whole spacetime as on the axis. Instead one has to deal with integral equations or to look 
directly for a differential relation. 

To avoid the differentiation of theta functions with respect to a branch point of the 
Riemann surface, we use the algebraic formulation of the hyperelliptic solutions (^) and 
(j39l). From the latter it can also be seen how one could get a relation independent of G 
without differentiation: one can consider the equations Q3"8|) and (KDJ) as integral equations 



for G. In principle one could try to eliminate G and X from these equations and (43). 
We will not investigate this approach but try to establish a differential relation. To this 
end it proves helpful to define the symmetric (in the K n ) functions Si via 

f[(K - Kt) =: K 9 - S r i^ rl + - + So (60) 
i=i 

i.e. So = K 1 K 2 . . . K g , . . . , S g _i = Ki + . . . + K g . The equations fl4"3| ) are bilinear in the 
real and imaginary parts of the Si which are denoted by Ri and respectively. With this 
notation we get 
Theorem 6.1: 

The Xi and the Ernst potential f are subject to the system of differential equations 
= (R -P R 1 + ... + P°(-iy)x z - 1 -Q 2 (P ) 

-~(1 - a; 2 )(ln/i% (l -P I 1 + ... + (61) 



2 



and for g > 1 



((-i) j+1 R j+1 + . . . + /;;' 1 •) - . • • • • • •''/''/ 

- x 2 )(\nff)M-iy +1 I J+1 + ...- Pt j -%-i)- (62) 

Proof: 

Differentiating (|43|) with respect to z and eliminating the derivatives of the K i z via 
the Picard-Fuchs relations (0), we end up with a linear system of equations for the 
derivatives of the and x which can be solved in standard manner. The Vandemonde- 
type determinants can be expressed via the symmetric functions. For x z one gets (|6l"D . 
The equations for the bilinear in the symmetric functions. They can be combined 

with to ©. 

Remark: If one can solve ( pE5D for the Ki, the equations QB"T| ) and ( p2] ) will be a non-linear 
differential system in z (and z which follows from the reality properties) for the Xj, x and 
/ which only contains the branch points of the Riemann surface as parameters. 
For the metric function a, we get with (|3"7|) 
Theorem 6.2: 

The metric function a is related to the functions Xi and Si via 

Z ~ " Vl " (63) 
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-h + (64) 



for g > 1 and 

i x(ai - C) 
1 -x 2 

forg = l. 
Proof: 

To express the function Z via the divisor X, we define the divisor T = T\ + . . . + T g as the 
solution of the Jacobi inversion problem uj(T) = u>(X) + w(-P) where P is in the vicinity 
of oo~ (only terms of first order in the local parameter near oo~ are needed). Using the 
formula for root functions (p3|), we get for the quantity Z in 



Z^^-lnfr^-^. (65) 

Applying Abel's theorem to the definition of T and expanding in the local parameter near 

oo~, we end up with (^) for general g > 1 and with ([54]) for g = 1. 

Remark: 

1. For g > 1 equation ( |6~3"|) can be used to replace the relation for x g _2 )Z in (|^) since the 
latter is identically fulfilled with (|63|) and (|TT|). 

2. An interesting limiting case is G ~ 1 where / ~ 1, i.e. the limit where the solution is 
close to Minkowski spacetime. By the definition fl3~8|), the divisor X is in this case approx- 
imately equal to D. Thus the symmetric functions in ( |62"D and (|6l|) can be considered as 
being constant and given by the branch points Ei. Relation ( |63D implies that the quantity 
Z is approximately equal to J 9 _2 in this limit, i.e. it is mainly equal to the constant ao 
in lowest order. Since the differential system flB"2] ) and ([H]) is linear in this limit, it is 
straight forward to establish two real differential equations of order g for the real and the 
imaginary part of the Ernst potential. In principle this works also in the non-linear case, 
where sign ambiguities in the solution of (|38D can be fixed by the Minkowskian limit. 
To illustrate the above equations we will first consider the elliptic case. This is the only 
case where one can establish an algebraic relation between Z and b independent of G. 
Equations ([43]) lead to 

(1 — x 2 )Rq = a,\ — (x 2 
(l-x 2 )S S = E 1 E 1 -P P x 2 . (66) 

Formula ( |64D takes with (^) (the sign of Iq is fixed by the condition that Iq = —f3\ for 
x = 0) the form 

(1 - x 2 )Z = ix(ax - C) + ^(1 - x 2 )(Pi - P 2 x 2 ) - x 2 ( ai - C) 2 . (67) 

This relation holds in the whole spacetime for all elliptic potentials, i.e. for all possible 
choices of G in (^). This implies that one can only solve boundary value problems on 
elliptic surfaces where the boundary data at some given contour T z satisfy condition floTI). 
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In the case genus 2, we get for 



(1 — x 2 )R\ = a.\ + a 2 — (x z + xx , 

{l-x 2 )(Rl + I 2 + 2R ) = (a 1 + a 2 ) 2 + 2a 1 a 2 + (3 2 + (3 2 -x 2 -x 2 (p 2 + ( 2 )+4:(xxo, 

(1 - x 2 )(R 1 R + Iil ) = aia 2 (a! + a 2 ) + ai(3 2 + a 2 /3 2 - Cxi + (p 2 + ( 2 )xx , 

(l-x 2 )(^ + / 2 ) = (a? + #)(af + #)-(^ + C 2 )4 (68) 



The aim is to determine the Si and Xq from (pq ) and 

(1 -x 2 )(Z + h) = ix -(ix, (69) 
and to eliminate these quantities in 

(R - P Ri + P^)x z = l -{x, + P x) + ~(1 - x 2 )(ln//),(/ - P h) (70) 

which follows from (|6l|). 

Remark: Boundary value problems 

Since the above relations will hold in the whole spacetime, it is possible to extend them 
to an arbitrary smooth boundary T z , where the Ernst potential may be singular (a jump 
discontinuity) and where one wants to prescribe boundary data (combinations of /, f z ). If 
these data are of sufficient differentiability (at least C 9,a (T z )), we can check the solvability 
of the problem on a given surface with the above formulas. The conditions on the differ- 
entiability of the boundary data can be relaxed by working directly with the equations 
(j38[) and (|39|) which can be considered as integral equations for InG. The latter is not 
very convenient if one wants to construct explicit solutions, but it makes it possible to 
treat boundary value problems where the boundary data are Holder continuous. We will 
only work with the differential relations and consider merely the derivatives tangential to 



T z in (|62"D to establish the desired differential relations between a, b and U. One ends 
up with two differential equations which involve only U, b and derivatives. The aim is to 
construct the spacetime which corresponds to the prescribed boundary data from these 
relations. To this end one has to integrate the differential relations using the boundary 
conditions. Integrating one of these equations, one gets g real integration constants which 
cannot be freely chosen since they arise from applying the tangential derivatives in (|62|) . 
Thus they have to be fixed in a way that the integrals on the right-hand side of (|38| ) 
are in fact the 6-periods of the second integral on the right-hand side of (|38|) and that 
(P9]) holds. The second differential equation arises from the use of normal derivatives of 
the Ernst potential in fl6~l"|). To satisfy the 6-period condition (0), one has to fix a free 
function in the integrated form of the corresponding differential equation. Thus one has 
to complement the two differential equations following from (|6T|) with an integral equa- 
tion which is obtained by eliminating G from e.g. u\ and u 2 in (^). For given boundary 
data, the system following from (^) may in principle be integrated to give e 2U and b 
in dependence of the boundary data. Then the in general non-linear integral equation 
will establish whether the boundary data are compatible with the considered Riemann 
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surface. This is typically a rather tedious procedure. There is however a class of prob- 
lems where it is unnecessary to use this integral equation. In case that the differential 
equations hold for an arbitrary function e 2U , the integral equation will only be used to 
determine this metric function, but the boundary value problem will be always solvable 
(locally). This offers a constructive approach to solve boundary value problems without 
having to consider non-linear integral equations. 



7 Counter-rotating disks of genus 2 

Since it is not very instructive to establish the differential relations for genus 2 in the 
general case, we will concentrate in this section on the form these equations take in the 
equatorially symmetric case for counter-rotating dust disks. In this case, the solutions are 
parametrized by E\ = a + i(3. We will always assume in the following that the boundary 
data are at least C 2 (T Z ) (the normal derivatives of the metric functions have a jump at 
the disk, but the tangential derivatives are supposed to exist up to at least second order). 
Putting s = be~ 2U and y = e 2U , we get for (|70f) for ( = 0, p < 1 



ix = (R Q -p 2 -sh) h ^-p{R 1 -sh) h ^-(s{R Q -p 2 ) + I )^ + p{sR l + h) y -^, 

V I y y V ) y y 

ps = (R Q -p 2 -sI ) b ^ + p(R 1 -sI 1 ) b ^-(s(R -p 2 )+I )^-p(sR 1 + I 1 )^, (71) 

V ) y y V ) y y 

where the Si and ixo are taken from (|68|) and (^). Since counter- rotating dust disks are 
subject to the boundary conditions (|23|), we can replace the normal derivatives in (|TT| ) 
via ([23|) which leads to a differential system where only tangential derivatives at the disk 
occur. With (BBT) and (p§|) we get 



• ^(7 ■ ^ R o-P 2 I R 2 + P 2 + 5y 2 Rp - p 2 \ f, y p 7 b p \ 

Z\ (Ro-P 2 R 2 + p 2 + Sy 2 \ ( v P , / 7 , • A X 



P s[l--) = l^-^^L- U S Z^ { -Z^ 1.(73) 



With h = ix /(l - x 2 ) - Z and 



\XqZ — a — ^ Z 2 — p 



2 



*>= ' 2^" (74) 

the function ixo follows from 

p2 , (Ro-P 2 ) 2 _ a 2 + P 2 -p 2 x 2 

Ko+ I 2 (l-x 2 ) 2 ~ 1-x 2 ' 1 j 

i.e. an algebraic equation of fourth order for ixo which can be uniquely solved by respect- 
ing the Minkowskian limit. Thus ( [72] ) and ( [73]) are in fact a differential system which 
determines b and y in dependence of the angular velocity Q. 
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7.1 Newtonian limit 



For illustration we will first study the Newtonian limit of the equations ([73]) (where 
counter-rotation does not play a role). This means we are looking for dust disks with 
an angular velocity of the form Q = ujq(p) where \q(p)\ < 1 for p < 1, and where the 
dimensionless constant uj « 1. Since we have put the radius p of the disk equal to 
1, uj — oopo is the upper limit for the velocity in the disk. The condition uj << 1 just 
means that the maximal velocity in the disk is much smaller than the velocity of light 
which is equal to 1 in the used units. An expansion in uj is thus equivalent to a standard 
post-Newtonian expansion. Of course there may be dust disks of genus 2 which do not 
have such a limit, but we will study in the following which constraints are imposed by the 
Riemann surface on the Newtonian limit of the disks where such a limit exists. 
The invariance of the metric fllQ|) under the transformation t — > — t and Q — > —Q implies 
that U is an even function in uj whereas b has to be odd. Since we have chosen an 
asymptotically non-rotating frame, we can make the ansatz y — 1 + uj 2 yi + • • ., b = 
uj 3 b 3 + . . ., and a = uj 3 a 3 + . . .. The boundary conditions (|23| ) imply in lowest order 
V2,p = 2g 2 p, the well-known Newtonian limit. Since (fP2|) reduces to the Laplace equation 
for y-i in order uj 2 , we can use the methods of section (|2|) to construct the corresponding 
solution. In order uj 3 , the boundary conditions ( |23"D lead to 

h, P = 2pqy2,o (76) 

whereas equation ( |T2"D leads to the Laplace equation for 63. Again we can use the methods 
of section |2|, but this time we have to construct a solution which is odd in £ because of the 
equatorial symmetry. In principle one can extend this perturbative approach to higher 
order, where the field equations (|12"D lead to Poisson equations with terms of lower order 
acting as source terms, and where the boundary conditions can also be obtained iteratively 
from (|23]). 

With this notation we get 
Theorem 7.1: 

Dust disks of genus 2 which have a Newtonian limit, i.e. a limit in which Q = ujq(p) 
where \q(p)\ < 1 for p < 1, are either rigidly rotating (q = 1) or q is a solution to the 
integro- differential equation 

b 3 =((R -p 2 )2q-K)y u (77) 

where in the first case I^/Rq = 2uj and in the second I\ = kuj with R® and k being uj 

independent constants. 

Proof: 

Since the right hand side of (l38f) vanishes, we have Ki = E{ for u> — > 0, and thus ao = I\ up 
to at least order uj 3 . Keeping only terms in lowest order and denoting the corresponding 
terms of the symmetric functions by S®, we obtain for (f73|) 

uj% = y u (2q{Rl - p 2 )uj 3 - uj 2 lf) . (78) 
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The second equation ( |T5D involves 63^ and is thus of higher order. If ( [75] ) holds, this 
equation will be automatically fulfilled. 

The ^-dependence in (|78|) implies that R®, J° and thus the branch points must depend 
on uj. Since ?/2,c is proportional to the density in the Newtonian case, it must not vanish 
identically. The possible cases following from equation (|78|) are constant Q or ([771). Using 
(^) and (|3|), one can express directly via Q which leads to 



IT JO p + p' 



with = 2y / p/7/(p + p'). Thus (|77|) is in fact an integro-differential equation for q. This 
completes the proof. 

7.2 Explicit solution for constant angular velocity and constant 
relative density 



The simplifications of the Newtonian equation (|78|) for constant f2 give rise to the hope 
that a generalization of rigid rotation to the relativistic case might be possible which we 
will check in the following. Constant j/Q makes it in fact possible to avoid the solution 
of a differential equation and leads thus to the simplest example. We restrict ourselves to 
the case of constant relative density, 7 = const. The structure of equation (|73| ) suggests 
that it is sensible to choose the constant a as a = —7/^ since in this case Z = R. 
This is the only freedom in the choice of the parameters a and (3 on the Riemann surface 
one has for g = 2 since one of the parameters will be fixed as in the Newtonian case by 
the condition that the disk has to be regular at its rim. The second parameter will be 
determined as an integration constant of the Picard-Fuchs system. 
We get 

Theorem 7.2: 

The boundary conditions ( \2$( ) and ( \2Q ) for the counter-rotating dust disk with constant 



n and constant 7 are satisfied by an Ernst potential of the form on a hyperelliptic 
Riemann surface of genus 2 with the branch points specified by 




« = -l + 2» /5 = V^ + (5 "T- (80) 
The parameter 5 varies between 5 = (only one component) and 5 = 5 S , 

Ss = 2 ^1 + \f+^) > (81) 
the static limit. The function G is given by 



G(r) = ^ • (82) 

^(r 2 - a) 2 + (3 2 - (r 2 + 1) 
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This is the result which was announced in [113 . 
Proof: 

The proof of the theorem is performed in several steps. 

1. Since the second factor on the right-hand side in flT3| ) must not vanish in the Newtonian 
limit, we find that for Z = R 



R -p 2 Z 2 + p 2 + 5y' 



h 2Z 
With this relation it is possible to solve ([74]) and (|69|) 



Z(p 2 + 2a-5y 2 (l-x 2 )) 
lx o = tto To T-o (84) 




^(1 - X 2 ) = —[i-S/j] + S[n -) ,\ [ A l(--nS) +<P.P. 



One may easily check that equation (|72|) is identically fulfilled with these settings. Thus 
the two differential equations flT2"| ) and (|73|) are satisfied for an unspecified y which implies 
that the boundary value problem for the rigidly rotating dust disk can be solved on a 
Riemann surface of genus 2 (the remaining integral equation which we will discuss below 
determines then y). 

2. To establish the integral equations which determine the function G and the metric 



potential e , we use equations (p8|). Since we have expressed above the Ki as a function 



of e alone, the left-hand sides of (BSI) are known in dependence of e . It proves helpful 



to make explicit use of the equatorial symmetry at the disk. By construction the Riemann 
surface £ is for ( = invariant under the involution K — > —K. This implies that the 
theta functions factorize and can be expressed via theta functions on the covered surface 
Ei given by p\(r) = r(r + p 2 )((r — a) 2 + (3 2 ) and the Prym variety S 2 (which is here also 
a Riemann surface) given by p^ij) = (r + p 2 )((r — a) 2 + (3 2 ) (see |]26], for details). On 
these surfaces we define the divisors V and W respectively via 

1 r-p 2 \nG(JT)dT [V dr 1 r^\nG(Jr)dr f w dr 
u v = — r~, =: / — , u w = — — =: / — . (85) 

17T JO / i l( r ) JO A 4 ! 17r J-P 2 /M r ) J°° ^2 



For the Ernst potential we get 




(86) 



where I v = 

3. Using Abel's theorem and (|38|), we can express V and W by the divisor X which leads 
to 2 

y = P X ° ($,7) 

ZHl-x 2 ) -2Zix n y ' 
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and 

1 

x 2 

4. Since V and I v vanish for p = 0, we can use (|86| ) for p = to determine the integration 



W + p 2 = -(Z 2 (l-x 2 )-2Zix -x 2 ) . (88) 



constant of the Picard-Fuchs system. We get with (|8 

1 „ 5 2 



{3 2 = --5a + -. (89) 



5. Since V in fl8T|) is with fl84] ) a rational function of p, a and /? and does not depend 
on the metric function e 2U , we can use the first equation in (p5|) to determine G as the 
solution of an Abelian integral which is obviously linear. With G determined in this way, 
the second equation in (|H) can then be used to calculate e 2U at the disk which leads to 
elliptic theta functions (see also |30|]). (In the general case, one would have to eliminate 
e 2U in the relations for u v and u w to end up with a non-linear integral equation for G.) 
The integral equation following from fl85|), 



v dr 1 r" 2 InG , 

dr (90) 



o Pi{t) m Jo pi[r 

is an Abelian equation and can be solved in standard manner by integrating both sides 
of the equation with a factor — r from to r where r = — p 2 . With (^7]) we get for 

what is essentially an integral over a rational function 



G(K) = y = 2 -. (91) 

'{K-a) 2 + f3 2 -(K-a + l) 

6. The condition G(—l) = 1 excludes ring singularities at the rim of the disk and leads 
to a continuous potential and density there. It determines the last degree of freedom in 



(0) to 

a = -l + 5 -- (92) 

7. The static limit of the counter-rotating disks is reached for (3 = 0, i.e. the value S s . 

This completes the proof. 

Remark: 

1. It is interesting to note that there are algebraic relations between a, b and e 2U though 
they are expressed via theta functions, i.e. transcendental functions, also at the disk. 

2. It is an interesting question whether there exist disks with non-constant 7/fi or S for 
genus 2 in the vicinity of the above class of solutions. Whereas this is rather straight 
forward for a non-constant S if 7/fi are constant, it is less obvious if the latter does not 
hold. This means that one looks for given S for solutions with 

1 = C + ep(p) (93) 
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where Co is a constant, \p\ < 1 is a function of p, and where e << 1 is a small dimensionless 
parameter. We can assume that p is not identically constant since this would only lead to 
a reparametrisation of the above solution. To check if there are solutions for small enough 
e, one has to redo the steps in the proof of theorem 7.2 in first order of e by expanding all 
quantities in the form y = y+ey+.... Doing this one recognizes that equation (|72|) becomes 
a linear first order differential equation for p of the form p p + F(p)p = where F is given 
by the solution for rigid rotation. For a solution p to this equation, the remaining steps 
can be performed as above. It seems possible to use the theorem on implicit functions 
to prove the existence of solutions for genus 2 in the vicinity of constant 7/O, but this is 
beyond the scope of this article. 



7.3 Global regularity 

In Theorem 7.2 it was shown that one can identify an Ernst potential on a genus 2 surface 
which takes the required boundary data at the disk. One has to notice however that this 
is only a local statement which does not ensure one has found the desired global solution 



which has to be regular in the whole spacetime except at the disk. It was shown in [29], £(| 
that this is the case if Q(u(oo~) + u) 7^ 0. In the Newtonian theory (see section (0)), 
the boundary value problem could be treated at the disk alone because of the regularity 
properties of the Poisson integral. Thus one knows that the above condition will hold in 
the Newtonian limit of the hyperelliptic solutions if the latter exists. For physical reasons 
it is however clear that this will not be the case for arbitrary values of the physical 
parameters: if more and more energy is concentrated in a region of spacetime, a black- 
hole is expected to form (see e.g. the hoop conjecture @]). The black-hole limit will be a 
stability limit for the above disk solutions. Thus one expects that additional singularities 
will occur in the spacetime if one goes beyond the black-hole limit. The task is to find 
the range of the physical parameters, here A and 5, where the solution is regular except 
at the disk. 
We can state the 
Theorem 7.3: 

Let E' be the Riemann surface given by p' 2 = (K 2 - E)(K 2 - E) and let a prime denote 
that the primed quantity is defined on E'. Let A c (<5) be the smallest positive value A for 
which e'(u') = 0. Then Q{uj{oo~) +u) ^ for all p, ( andO < A < X c (5) andO <S < 5 S . 
This defines the range of the physical parameters where the Ernst potential of Theorem 



7.2 is regular in the whole spacetime except at the disk. Since it was shown in |29|, [3C 
that O'(u') = defines the limit in which the solution can be interpreted as the extreme 
Kerr solution, the disk solution is regular up to the black-hole limit if this limit is reached. 
Proof: 

1. Using the divisor X of (|3^) and the vanishing condition for the Riemann theta function, 
we find that 0(u;(oo~) +u) = is equivalent to the condition that oo + is in X. The reality 
of the Ui implies that X = oo + + (— iz). Equation (^) thus leads to 

°° + dr f~ iz dr 1 f In Gdr 

— + 73 7 / = 0, 

Ei p Je 2 p 2tti Jr p 
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rdr f-' lz rdr 1 r InGrdr , , 

+ / —-t^ — = 0, (94) 



e x p Je 2 p 2-iri Jr p 

where = denotes equality up to periods. The reality and the symmetry with respect to £ 
of the above expressions limits the possible choices of the periods. It is straight forward 
to show that 6(c<j(oo~) + u) = if and only if the functions F defined by 

Fx := I -+I ~--nx\2i ^ + 2 /-+/ - + 




F 2 := I —+, ^--n 2 \2i T *L + 2 I T *L + 



fel /l Jb 2 A'' Jdl p >' a 2 A* 

f ^ + 2/ Z* 

bi /(X ^2 A* 




(95) 



with the cut system of Fig. 1 and with 711,2 £ Z vanish for the same values of p, £, A, 5. 
The functions F are both real, F\ is even in ( whereas F 2 is odd. Thus F 2 is identically 
zero in the equatorial plane outside the disk. 

2. In the Newtonian limit A ~ 0, the above expressions take in leading order of A the 
form 

Fx = A ((-8rn + l) Cl (p, C) In A - d x {p, ()\) , (96) 

and 

F 2 = VX ((-8n 2 + l)c 2 (p, C) In A - d a (p, C)A§) , (97) 



where we have used the same approach as in the calculation of the axis potential in 
(see |30 1 and references given therein); the functions cx, dx are non- negative whereas c 2 /d 2 
is positive in C/{C = 0}. Thus the F are zero for A = which is Minkowski spacetime 
f — 1, but they are not simultaneously zero for small enough A, i.e. / is regular in the 
Newtonian regime in accordance with the regularity properties of the Poisson integral. 
The Fi may vanish however at some value A s for given p, ( and 5. Since we are looking 
for zeros of the F in the vicinity of the Newtonian regime, we may put nx,2 = 1 here. 

3. Let Q be the open domain C/{( = 0, p < 1 Vp = 0}. It is straight forward to check that 
the F are a solution to the Laplace equation AFj = with A = 4 (d z2 + 2 (^+ g ) + ^z)) 
for z, z G Q. Thus by the maximum principle the F, do not have an extremum in Q. 

4. At the axis for Q > 0, the Uj are finite whereas the Fj diverge proportional to — In p for 
all A, 5. Thus / is always regular at the axis. 

5. Relation (|6l|) at the disk can be written in the form (y + A) 2 + b 2 = B 2 where A and 
B are finite real quantities. Thus the Ernst potential is always regular at the disk. Due 
to symmetry reasons F 2 = u 2 which is non-zero except at the rim of the disk. For Fx one 
gets at the disk 

r°° + dr r E dr r B dr . . 

Fx = / — - + / — - + / — - - u v . 98 
J- P 2 n x [j) Jo pxij) Jo px(r) 

With ( p0|) one can see that Fx is always positive at the disk. 

6. Since Fx is strictly positive on the axis and the disk and a solution to the Laplace 
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equation in Q, it is positive in C if it is positive at infinity. F\ is regular for \z\ — > oo and 
can be expanded as F\ = Fn/\z\ + o(l/\z\) where Fn can be expressed via quantities on 
E'. We get 

F * & ^ T ^ t ^ n( ^ r 
11 ~ 2 /&; /7 ~ 2vri 7_i // ' ^ 

The quantity Fn = iff Q'(u f ) = 0. The condition Fn > is thus equivalent to the 
condition that A < A c (<5) where X C (S) is the first positive zero of O'(tt'). This completes 
the proof. 
Remark: 

1. In the second part of the paper we will show that the ultrarelativistic limit (vanishing 
central redshift) in the case of a disk with one component is given by O'(w') = for a 
finite value of A. In the presence of counter-rotating matter, this limit is however not 
reached, the central redshift diverges for A = oo and Q'(u') ^ 0. This supports the 
intuitive reasoning that counter-rotation makes the solution more static, i.e. it behaves 
more like a solution of the Laplace equation with the regularity properties of the Poisson 
integral. 

2. Since F%(p, 0) = for p > 1, the reasoning in 6. of the above proof shows that there will 
be a zero of 6(u;(oo~) + u) and thus a pole of the Ernst potential in the equatorial plane 
for A > X C (S) if the theta function in the numerator does not vanish at the same point. In 
the equatorial plane the Ernst potential can be expressed via elliptic theta functions (see 

30fl ) which have first order zeros. Thus Fn will be negative for A > A c in the vicinity of 



A c , and consequently the same holds for F\ in the equatorial plane at some value p > 1. It 
will be shown in the third article that the spacetime has a singular ring in the equatorial 
plane in this case. The disk is however still regular and the imposed boundary conditions 
are still satisfied. This provides a striking example that one cannot treat boundary value 
problems locally at the disk alone in the relativistic case. Instead one has to identify the 
range of the physical parameters where the solution is regular except at the disk. 



8 Conclusion 

We have shown in this paper how methods from algebraic geometry can be successfully 
applied to construct explicit solutions for boundary value problems to the Ernst equation. 
We have argued that there will be differentially rotating dust disks for genus 2 of the 
underlying Riemann surface in addition to the one we could identify explicitly. To prove 



existence theorems for solutions to boundary value problems, the methods of [II], f|2 
seem to be better suited since the hyperelliptic techniques are limited to finite genus of the 
Riemann surface. Moreover the used techniques at the boundary have to be complemented 
by a proof of global regularity. The finite genus of the Riemann surface also restricts the 



usefulness in the numerical treatment of boundary value problems. The methods of [12] 



and jT3j] are not limited in a similar way and have proven to be highly efficient. Thus 
the real strength of the approach we have presented here is the possibility to construct 
explicit solutions whose physical features can then be discussed in analytic dependence 
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on the physical parameters up to the ultrarelativistic limit. Whether this approach can 
be generalized to more sophisticated matter models or whether the equations can still be 
handled for higher genus will be the subject of further research. 
Acknowledgement 

I thank J. Frauendiener, R. Kerner, D. Korotkin, H. Pfister, O. Richter and U. Schaudt 
for helpful remarks and hints. The work was supported by the DFG and the Marie-Curie 
program of the European Community. 

References 

[1] J. B. Hartle and D. H. Sharp, Astrophys. J., 147, 317 (1967). 

[2] L. Lindblom, Astrophys. J., 208, 873 (1976). 

[3] H. D. Wahlquist, Phys. Rev., 172, 1291 (1968). 

[4] M. Bradley, G. Fodor, M. Marklund and Z. Perjes, Class. Quant. Grav., 17, 351 
(2000). 

[5] F. J. Ernst, Phys. Rev, 167, 1175; Phys. Rev 168, 1415 (1968). 

[6] D. Maison, Phys. Rev. Lett, 41, 521 (1978). 

[7] V. A. Belinski and V. E. Zakharov, Sov. Phys.-JETP, 48, 985 (1978). 

[8] G. Neugebauer, J. Phys. A: Math. Gen., 12, L67 (1979). 

[9] C. Struck, Physics Reports, 321, 1 (1999). 

[10] J. Bicak, D. Lynden-Bell and J. Katz, Phys. Rev. D, 47, 4334 (1993); MNRAS, 265, 
126 (1993). 

[11] J. Bicak and T. Ledvinka, Phys. Rev. Lett, 71, 1669 (1993). 

[12] T. Nozawa, N. Stergioulas, E. Gourgoulhon, Y. Eriguchi, Astron. Astrophys. Suppl. 
Ser., 132, 431 (1998). 

[13] A. Lanza in Approaches to Numerical Relativity, ed. by R. D'Inverno (Cambridge: 
Cambridge University Press) 281 (1992). 

[14] J. M. Bardeen and R. V. Wagoner, Ap. J., 167, 359 (1971). 

[15] J. Binney and S. Tremaine, Galactic Dynamics (Princeton Univ. Press, Princeton, 
1987). 

[16] C. Klein and O. Richter, J. Geom. Phys., 24, 53, (1997). 
[17] C. Klein and O. Richter, Phys. Rev. D, 57, 857 (1998). 



30 



[18] D. A. Korotkin, Theor. Math. Phys. 77, 1018 (1989), Commun. Math. Phys., 137, 
383 (1991), Class. Quant. Grav., 10, 2587 (1993). 

[19] C. Klein and O. Richter, Phys. Rev. Lett, 83, 2884 (1999). 

[20] T. Morgan and L. Morgan Phys. Rev., 183, 1097 (1969); Errata: 188, 2544 (1969). 

[21] G. Neugebauer and R. Meinel, Astrophys. J., 414, L97 (1993), Phys. Rev. Lett., 
73, 2166 (1994), Phys. Rev. Lett, 75, 3046 (1995). 

[22] D. Kramer, H. Stephani, E. Herlt and M. MacCallum, Exact Solutions of Einstein's 
Field Equations, Cambridge: CUP (1980). 

[23] C. Pichon and D. Lynden-Bell, MNRAS, 280, 1007 (1996). 

[24] W. Israel, Nuovo Cimento 44B 1 (1966); Errata: Nuovo Cimento, 48B, 463 
(1967). 

[25] H. M. Farkas and I. Kra, Riemann surfaces, Graduate Texts in Mathematics, Vol. 71, 
Berlin: Springer (1993). 

[26] E.D. Belokolos, A.I. Bobenko, V.Z. Enolskii, A.R. Its and V.B. Matveev, Algebro- 
Geometric Approach to Nonlinear Integrable Equations, Berlin: Springer, (1994). 

[27] D. Korotkin and V. Matveev, Fund. Anal. Appl, 34 1 (2000). 

[28] D. Korotkin and V. Matveev, Lett. Math. Phys., 49, 145 (1999). 

[29] C. Klein and O. Richter, Phys. Rev. Lett, 79, 565 (1997). 

[30] C. Klein and O. Richter, Phys. Rev. D, 58, CID 124018 (1998). 

[31] R. Meinel and G. Neugebauer, Phys. Lett. A 210, 160 (1996). 

[32] D. A. Korotkin, Phys. Lett. A 229, 195 (1997). 

[33] M. Ansorg and R. Meinel, kr-qc/9910045l (1999). 

[34] R. Geroch J. Math. Phys. 11, 2580 (1970). 

[35] R. Hansen J. Math. Phys. 15, 46 (1974). 

[36] G. Fodor, C. Hoenselaers and Z. Perjes, J. Math. Phys., 30, 2252 (1989). 

[37] I. Hauser and F. Ernst, J. Math. Phys., 21, 1126 (1980). 

[38] P. Kordas, Class. Quantum Grav., 12, 2037, (1995). 

[39] R. Meinel and G. Neugebauer, Class. Quantum Grav., 12, 2045, (1995). 



31 



[40] K. Thorne in Magic Without Magic: John Archibald Wheeler, ed. by J. Klauder 
(W.H. Freeman: San Francisco), 231 (1972). 

[41] U. Schaudt and H. Pfister, Phys. Rev. Lett, 77, 3284 (1996) 

[42] U. Schaudt, Comm. Math. Phys., 190, 509, (1998) 



32 



